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Spatial variations in the distribution of galaxy luminosities, estimated from redshifts as distance 
proxies, are correlated with the peculiar velocity field. Comparing these variations with the peculiar 
velocities inferred from galaxy redshift surveys is a powerful test of gravity and dark energy theories 
on cosmological scales. Using ~ 2 x 10® galaxies from the SDSS Data Release 7, we perform this 
test in the framework of gravitational instability to estimate the normalized growth rate of density 
perturbations fag. = 0.37 ±0.13 at ^ ~ 0.1, which is in agreement with the ACDM scenario. This 
unique measurement is complementary to those obtained with more traditional methods, including 
clustering analysis. The estimated accuracy at 2 ~ 0.1 is competitive with other methods when 
applied to similar datasets. 
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Introduction .—Unraveling the origin of cosmic acceler¬ 
ation remains one of the biggest challenges in fundamen¬ 
tal physics. Lacking a natural explanation within the 
standard paradigms of cosmology and particle physics, 
many theoretical models have been proposed, ranging 
from the inclusion of new scalar fields to genuine modi¬ 
fications of general relativity [l|, 0 . Models that provide 
the same expansion history of the universe can generally 
lead to a very different evolution of density fluctuations. 
In the linear regime, the latter is fully captured by the 
growth rate. 
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where SI is the matter density parameter, a is the cosmic 
scale factor, and D{a) is the growing mode of density per¬ 
turbations. The standard cold dark matter model with a 
cosmological constant (ACDM) is well characterized by 
a growth index 7 = dlog//dlogSI ~ 6/11 Q, but other 
frameworks typically yield different values and may even 
exhibit a scale-dependent behavior [e.g., ii- Sensitive 
to both expansion and perturbations, constraining / thus 
represents a key observational test to disentangle compet¬ 
ing world models, which justifies the increasing number 
of projects aimed at its measurement. 

Among the most robust probes of / is the analysis of 
redshift-space distortions (RSDs) p-lsj, i.e. the appar¬ 
ent anisotropic clustering of extragalactic objects such 
as galaxies |9l4l2l| and quasars |13l4l5l| in redshift space, 
with an accuracy mainly limited by the uncertainty of 
modeling nonlinearities in gravitational dynamics. Meth¬ 
ods based on direct measurements of peculiar velocities 
[AS have also been successful, but the sparseness and 
small galaxy numbers in these catalogs restrict them to 
the very local universe at redshifts 2 ~ 0. Other sug¬ 
gested app roaches include weak leasing, cluster abun¬ 
dances [Si, and Lyman-alpha absorption [l^ . 


In this letter, we apply an independent method, based 
on luminosity modulations, to galaxies from the Sloan 
Digital Sky Survey Data Release 7 (SDSS DR7) [ 2 ^ 211, 
and obtain a new estimate of the growth rate at 2 0 . 1 . 

Altered by the line-of-sight components of their peculiar 
motions, v, galaxy redshifts generally differ from their 
actual distances |22l ] . Consequently, intrinsic magnitudes 
M inferred from the observed flux using redshifts appear 
brighter or dimmer than their true values, 
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where Dl is the luminosity distance, Zc denotes the cos¬ 
mological redshift due to the Hubble expansion, and the 
superscript (t) indicates the true quantities that would 
be measured if galaxies were correctly placed at Zc rather 
than 2 . For nearby galaxies and vjcz 1 (c is the speed 
of light), the contribution of velocities to the magnitude 
shift amounts to M — ^ 2.17vjcz. On large scales, 

the peculiar velocity field is spatially coherent and inde¬ 
pendently inferable from the observed galaxy distribution 
in space. Gravitational instability theory [ 2 ^ predicts 
the velocity field as a function of /3 = f/b, where b is 
the linear bias between galaxies and the total mass. This 
prediction allows a correction to M, and hence, bounds 
on /3 can be found by minimizing the scatter of corrected 
magnitudes with respect to a reference distribution 241. 

Data .—The analysis is based on the latest version of 
the NYU Value-Added Galaxy Catalog (NYU-VAGC) 
derived from the SDSS DR7 [261 ]. To minimize incom¬ 
pleteness and systematics, we focused on the subsample 
safe and considered the same dataset as in Ref. [l^. We 
then selected only galaxies in the range 0.06 < 2 < 0.12 
inside a rectangular patch specified in SDSS survey co¬ 
ordinates, —33° < 77 < 36° and —48° < A < 51.5°. 
This luminosity sample, which we used in the likeli¬ 
hood analysis below, is flux-limited in the r-band and 
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FIG. 1. Left .—Histograms of /ag-estimates (and Gaussian fits) derived with the luminosity approach from 128 mock catalogs, 
considering only multipoles I > 1 (solid lines) and I > 5 (dashed lines) in the velocity reconstructions. The distributions peak 
near the true value, (/o'8)true ~ 0.44, and deviate from a symmetric Gaussian mainly because of the velocity field’s nonlinear 
/3-dependence. Right .—Estimated Ax^ (and quadratic approximations) as a function of fag for real SDSS galaxies and velocity 
models with I > 1 (filled circles) and I > 5 (open circles). The results assume the measured power spectrum amplitude of 
L*-galaxies from Ref. [2^. 


contains around 2 x 10^ galaxies. To build peculiar ve¬ 
locity models, we also considered a second sample, tagged 
velocity, which corresponds to the maximum volume- 
limited sample within 0.05 < z < 0.13 trimmed to the 
range 0.06 < z < 0.12 and typically includes ~ 8 x 10“* 
galaxies for spatially flat cosmologies with 12 ss 0.3. To 
compute absolute magnitudes, we assumed a linear evo¬ 
lution model Q(z) = 1.6(z — 0.1) [13 and iT-corrections 
from the NYU-VAGC [H]. 

In addition, we used realistic mock SDSS catalogs gen¬ 
erated from the Millennium Simulation (^ . to cali¬ 
brate the velocity models and to assess the uncertainty on 
the /3-estimates. These mocks were customized to match 
data characteristics such as number counts, luminosity 
distribution, and sky coverage. 

Method .—A full account of the luminosity method can 
be found in Ref. 2^. Here we provide a brief summary 
of its key ingredients. Given a galaxy survey with mag¬ 
nitudes, spectroscopic redshifts, and angular positions 
on the sky, one traces the 3D galaxy distribution and re¬ 
constructs the liuear peculiar velocity field as a function 
of /3 [3l|. The correct value of /3 is then estimated by 
maximizing the probability of observing the data, 


Ptot = \{p {M,\z,,vim 
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ity field evaluated at the position of galaxy i, and redshift 
errors (typically about 10“"^ for SDSS galaxies, yielding 
sharp sample cuts in redshift space) are neglected 
Here 4>{M) is the galaxy luminosity function (LF) deter¬ 
mined from the full dataset, and the limiting magnitudes 
depend on u(/3) through the cosmological redshift 
Zc- The goal of this approach is to find the /3-value which 
minimizes the spread in the observed magnitudes. 

Our study does not depend on the Hubble constant and 
assumes a AGDM cosmology with fixed density parame¬ 
ters taken from Ref. 33| . Following the procedure of Ref. 
31j . the velocities were reconstructed in spherical har¬ 


where Ui(/3) denotes the radial part of the peculiar veloc¬ 


monics up to a multipole /max = 150 after smoothing the 
galaxy density field with a Gaussian kernel of 10h~^ Mpc 
radius. The monopole and dipole contributions cannot be 
reliably modeled for the dataset considered here. Exclud¬ 
ing them from the velocity model (/ > 1), we worked with 
a full sphere and fixed the boundary conditions by set¬ 
ting the density contrast outside the observed volume to 
zero. This particular approach can bias the low-/ modes 
of the velocity model, but the effect is easily estimated 
and removed with the help of the mock galaxy catalogs. 
Since the data covers a narrow z-range, we assumed no 
evolution of /3 and b in the analysis, and computed a set 
of model velocity fields varying /3 between 0 and 1 in steps 
of A/3 = 0.05. Dominated by large-scale structure, the 
cosmic velocity field is iusensitive to small-scale features 
such as those related to galaxy bias. We checked that 
our models aud results are robust to the precise choice of 
the smoothing scale. 
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Assigning v{j3) to the subsample luminosity, we re¬ 
moved galaxies near the sample edges (around 10%) to 
avoid artifacts due to incorrect boundary conditions. To 
maximize the probability Ptot j we used a spline-based 
LF estimator with a separation AM = 0.5 27, 34|. The 
partial sky coverage yields a statistical mixing between 
the velocity models and the LF parameters, which may 
result in biased constraints on /3. To deal with this is¬ 
sue, we set the LF to its estimate for a vanishing velocity 
field, evaluated Ptot for the different models of u(/3), and 
determined the maximum probability through interpola¬ 
tion. The corresponding bias on /3 was then inferred from 
the mock catalogs and is at the level of few percent. Al¬ 
though acceptable in view of the statistical uncertainties, 
we corrected for this bias in our analysis. 

Results. —The left panel of Fig. [T] shows the distribu¬ 
tion of growth rates estimated from 128 mock catalogs 
for velocity models with I > 1. We express our results 
in terms of fas = Paf^^, where is the measured 
amplitude of galaxy number counts in spheres of 8h~^ 
Mpc radius, related to the amplitude of mass fluctua¬ 
tions by as = af^^/b. This combination eliminates b in 
the comparison to other measurements and cosmological 
predictions [s^. In addition to the case I > 1, we also 
considered models with / > 5 to test the robustness to 
low multipoles which are most susceptible to the choice 
of boundary conditions. In both cases, the distributions 
peak around the true value, {fas)true ~ 0.44, showing an 
average spread of about 40%. The results deviate from a 
symmetric Gaussian distribution mainly because of the 
velocity field’s nonlinear dependence on (3. Shifting the 
multipole cutoff from I > 1 to I > 5 removes information 
and slightly increases the observed spread as expected. 

The analysis of the real data yielded f3 = 0.42 ± 0.14 
for 1 > 1, which, using the measured power spectrum 
amplitude of P*-galaxies from Ref. [2^, translates into 
fas = 0.37±0.13 (see right panel of Fig. [T]). The quoted 
errors were derived from the quadratic approximation of 
the log-likelihood around its maximum. Similarly, we es¬ 
timated f} = 0.63 ± 0.28 and fas = 0.56 ± 0.25 for I > 5. 
The corresponding relative errors of 35% and 45% are 
consistent with the mean scatter in Fig. [T]and the spread 
of errors computed for individual mocks. Both estimates 
are stable against removing high-1 modes (1 > 100) and 
agree well with the standard ACDM model which pre¬ 
dicts a value of fas ~ 0.42 for our fiducial cosmology. 
Other choices of cosmological parameters have no signifi¬ 
cant impact on these results. Setting 0 = 1, for instance, 
leads to changes of around 13% for / > 1 and 4% for I > 5. 

Performing a suite of basic tests 27|, we have also 
verified that our results are insensitive to the adopted 
AT-corrections and the modeling of luminosity evolution. 
Photometric uncertainties and a possible environmen¬ 
tal dependence of the LF play a minor role since the 
method, comparing luminosities to peculiar velocities ob¬ 
tained from the very same dataset, is remarkably robust 



FIG. 2. Comparison of growth rate measurements at low 
redshifts. Shown are constraints on fas from different red- 
shift surveys based on RSDs pH . I37l - l4(ij| , galaxy luminosities 
(filled square and circle) [IJ, and direct estimates of pecu¬ 
liar velocities (dashed, empty square) [l^ . All measurements 
assume a fixed cosmology and the error bars indicate the re¬ 
spective 68% confidence intervals. The shaded areas denote 
the 68% and 95% confidence limits inferred from the Planck 
data (TT-flowP-flensing) [^ . The square data points corre¬ 
spond to 2 = 0.02 and are offset for clarity. 


to these effects 

Conclusions .—Spatial modulations in the distribution 
of estimated galaxy luminosities trace the cosmic peculiar 
velocity field. As demonstrated by preliminary studies of 
very local samples m 0] , this can be combined with 
velocity reconstruction techniques to probe the linear 
growth factor with galaxy redshift surveys. The modula¬ 
tions in the luminosities of ~ 2 x 10® SDSS galaxies from 
the NYU-VAGC yield a value of fas at z ~ 0.1 which 
is in agreement with the ACDM cosmological model as 
dictated by the Planck data [s^. 

In Fig. [21 we present a compilation of recent measure¬ 
ments of the growth rate from different low-redshift sam¬ 
ples (z < 0.3) which were obtained from RSDs [HI, 
l40|, galaxy luminosities (filled square and circle) [24|, 
and direct estimates of peculiar velocities (dashed, empty 
square) [I^. The consistency between these measure¬ 
ments is striking in view of the different possible system¬ 
atic biases associated with the three independent meth¬ 
ods. At z 0.1, the error in the estimate of fas given 
here is competitive with that of other methods. The ac¬ 
curacy of our estimate (filled circle) matches that from 
the RSD analysis of SDSS main galaxies conducted in 
Ref. [i^ (empty triangle) which adopts a dataset with 
similar characteristics. Methods based on luminosities or 
distance indicators such as the Tully-Fisher relation M| 
are less sensitive to nonlinear corrections than the two- 
point statistics which enter the analysis of RSDs. Gom- 
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bini ng these approaches will result in superior control 

:il potential svstematics and allow improved bounds 


ove: 


on 


the growth rate. Further, this will help to resolve the ob- 
serv ed trend in RSD measurements of /gg which are con - 
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